Amplification of Fluctuations in a Spinor Bose Einstein Condensate 
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Dynamical instabilities due to spin-mixing collisions in a 87 Rb F = 1 spinor Bose-Einstein con- 
densate are used as an amplifier of quantum spin fluctuations. We demonstrate the spectrum of this 
amplifier to be tunable, in quantitative agreement with mean-field calculations. We quantify the 
microscopic spin fluctuations of the initially paramagnetic condensate by applying this amplifier and 
measuring the resulting macroscopic magnetization. The magnitude of these fluctuations is consis- 
tent with predictions of a beyond-mean-field theory. The spinor-condensate-based spin amplifier is 
thus shown to be nearly quantum-limited at a gain as high as 30 dB. 



1.2 



Accompanied by a precise theoretical framework and 
created in the lab in a highly controlled manner, ultracold 
atomic systems serve as a platform for studies of quan- 
tum dynamics and many-body quantum phases. Among 
these systems, gaseous spinor Bose Einstein condensates 
[E 0, 13,0, H> m which atoms may explore all sub-levels 
of a non-zero hyperfine spin F, provide a compelling op- 
portunity to access the static and dynamical properties 
of a magnetic superfluid [H, 0, H, [tj E3] • 

We previously identified a quantum phase transition 
in an F = 1 spinor Bose Einstein condensate between a 
paramagnetic and ferromagnetic phase 0]. This transi- 
tion is crossed as the quadratic Zeeman energy term, of 
the form qF 2 , is tuned through a critical value q = qo; 
here, F z is the longitudinal (z axis) projection of the di- 
mensionless vector spin operator F. Accompanying this 
phase transition is the onset of a dynamical instability 
in a condensate prepared in the paramagnetic ground 
state, with macroscopic occupation of the \m z = 0} mag- 
netic sublevel 11, 12, r|. This instability causes trans- 
verse spin perturbations to grow exponentially, produc- 
ing atoms into the \m z — ±1) sublevels. In contradic- 
tion with the mean-field prediction that the paramag- 
netic state should remain stationary because it lacks fluc- 
tuations by which to seed the instability, experiments 
revealed the spontaneous magnetization of such conden- 
sates after they were rapidly quenched across the phase 
transition. 

In this Letter, we investigate the spin fluctuations that 
become amplified by the spin-mixing instability. In par- 
ticular, we test whether these fluctuations correspond to 
quantum noise, i.e. to the zero-point fluctuations of quan- 
tized spin excitation modes that become unstable. For 
this, we use the spin-mixing instability as an amplifier, 
evolving microscopic quantum fluctuations into measur- 
able macroscopic magnetization patterns. We present 
two main results. First, we characterize the spin-mixing 
amplifier and demonstrate its spectrum to be tunable 
by varying the quadratic Zeeman shift. This spectrum 
compares well with a theoretical model that accounts for 
the inhomogeneous condensate density and for magnetic 
dipole interactions. Second, we measure precisely the 



transverse magnetization produced by this amplifier at 
various stages of amplification, up to a gain of 30 dB in 
the magnetization variance. This magnetization signal 
corresponds to the amplification of initial fluctuations 
with a variance slightly greater than that expected for 
zero-point fluctuations. 

Descriptions of the dynamics of initially paramagnetic 
spinor condensates [HI, [H, Hi , HI, 3 13] have focused 
on the effects of the quadratic Zeeman shift and of the 
spin-dependent contact interaction. The latter interac- 
tion has the mean-field energy density C2«(F) 2 , and, with 
c 2 = 4:Trh 2 Aa/3m < 0, favors a ferromagnetic state; here, 
Aa = (a2— do) where ajr tot is the s-wave scattering length 
for collisions between particles of total spin F to t [HQ and 
m is the atomic mass. Excitations of the uniform con- 
densate include both the scalar density excitations and 
also two polarizations of spin excitations with a disper- 
sion relation given as Eg (k) = (sk + q)(£k + q — qo), 
where Sk — h 2 k 2 /2m and qo — 2c2ti. For q > qo, spin 
excitations are gapped (E 2 > 0), and the paramagnetic 
condensate is stable. Below this critical value, the para- 
magnetic phase develops dynamical instabilities, defined 
by the condition E 2 < 0, that amplify transverse magne- 
tization. The dispersion relation defines the spectrum of 
this amplification, yielding a wavevector-dependent time 
constant for exponential growth of the power in the un- 
stable modes, r = h/2^/\E 2 \. 

The unstable regime is divided further into two regions. 
Near the critical point, reached by a "shallow" quench to 
<Zo/2 < q < qo, the fastest-growing instability occurs at 
zero wavevector, favoring the "light-cone" evolution of 
magnetization correlations at ever- longer range ll| . For 
a "deep" quench, with q < qo/2, the instabilities reach 
a maximum growth rate of 1/r = qo/h. The non-zero 
wavevector of this dominant instability sets the size of 
magnetization domains produced following the quench. 

Experimentally, we characterize the spectrum of this 
amplifier by seeding it with broadband noise and then 
measuring precisely the spectrum of its output. Sim- 
ilar to previous work we produce condensates of 
No = 2.0 x 10 6 87 Rb atoms, with a peak density of 
n = 2.6(1) x 10 14 cm -3 and a kinetic temperature of 



2 




~ 50 nK, trapped in a linearly polarized optical dipole 
trap characterized by trap frequencies (u) T ., lo,,, oj z ) = 
2ttx(39,440,4.2)s- 1 . Taking Aa = -1.4(3) a B with 
being the Bohr radius, the spin healing length £ s = 
(87rn|Aa|) -1 / 2 = 2.5 /im is larger than the condensate ra- 
dius r y — 1.6 (im along the imaging axis (y) . Thus, the 
condensate is effectively two-dimensional with respect to 
spin dynamics. For this sample, go — 2c2(n) = hx 15 Hz 
given the maximum y-axis column-averaged condensate 
density n. 

The quadratic Zeeman shift arises from the applica- 
tion of both static and modulated magnetic fields. A 
constant field of magnitude B, directed along the long 
axis of the condensate, leads to a quadratic shift of 
<lB/h = (70Hz/G 2 )£> 2 . In addition, a linearly polarized 
microwave field, detuned by S/2ir — ±35 kHz from the 
\F = 1, m z = 0) to \F = 2, m z = 0) hyperfine transition, 
induces a quadratic (AC) Zeeman shift of = —Ml 2 /AS 
where f2 is the Rabi frequency for the driving field. 

The condensate is prepared in the \m z = 0) state us- 
ing rf pulses followed by application of a 6 G/cm magnetic 
field gradient that expels atoms in the \m z = ±1) states 
from the optical trap 0]. This preparation takes place 
in a static 4 G field and with no microwave irradiation, 
setting q = qb + q^ > qa so that the paramagnetic con- 
densate is stable. Next, we increase the microwave field 
strength to a constant value, corresponding to a Rabi fre- 
quency in the range of 2ir x (0 - 1.5) kHz, to set q^. To 
initiate the instability, we ramp the magnetic field over 
5 ms to a value of B = 230 mG (giving qs/h = 7.6 Hz). 
During separate repetitions of the experiment (for dif- 
ferent values of the quadratic Zeeman shift at the 
end of the ramp was thus brought to final values qf/h 
between -2 and 16 Hz. 

Following the quench, the condensate spontaneously 
develops macroscopic transverse magnetization, saturat- 
ing within about 110 ms to a pattern of spin domains, 
textures, vortices and domain walls Q. Using a 2-ms- 
long sequence of phase-contrast images, we obtain a de- 
tailed map of the column-integrated magnetization M at 
a given time after the quench 3, 20| . The experiment 
is then repeated with a new sample. 

The observed transverse magnetization profiles of 
spinor condensates (Fig. QJ confirm the salient features 



predicted for the spin- mixing amplifier 27j. The vari- 
ation of the amplifier's spatial spectrum with qf is re- 
flected in the characteristic size Id of the observed spin 
domains, taken as the distance from the origin at which 
the magnetization correlation function, 
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acquires its first minimum; here gF^B is the atomic mag- 
netic moment. This characteristic size increases with in- 
creasing qf (Fig. 2). For qj/h > 9 Hz, the magnetization 
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FIG. 1: Transverse magnetization produced near the con- 
densate center after 87 ms (top) and 157 ms (bottom) of 
amplification at variable quadratic Zeeman shift qf. Mag- 
netization orientation is indicated by hue and amplitude by 
brightness (color wheel shown). The characteristic spin do- 
main size grows as qf increases. The reduced signal strength 
for qf/h>9 Hz reveals the lower temporal gain of the spin- 
mixing amplifier. 



features become long ranged and, thus, dominated by 
residual magnetic field inhomogeneities (< 2/iG). 

The data also confirm the distinction between deep and 
shallow quenches. The spatially averaged magnetization 
strength during the amplification, quantified by G(0) at 
t = 87 ms after the quench, is found to be constant for 
< qf/h < 6 Hz, reflecting that the temporal gain of the 
amplifier is uniform over < q < qo/2 [ill- For shallow 
quenches, with qf/h > 7 Hz, the measured magnetization 
decreases, reflecting a diminishing gain with increasing qf 
up to the transition point. 

While the above observations are consistent with the- 
oretical predictions, we note the unexpected outcome of 
quenches to negative values of qf. Such quenches revealed 
a diminished amplifier gain, resulting in a complete sup- 
pression of the growth of magnetization for qf/h< —1 
Hz. We are unable to account for this behavior. 

Having characterized the spin- mixing amplifier, let us 
consider the source of its input signal. For this, we de- 
velop a quantum field description of the spin-mixing in- 
stability ll|, working in the polar spin basis, where (\> n ^ 
is the annihilation operator for atoms of wavevector k 
in the zero-eigenvalue states of F ■ n. Treating a uni- 
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FIG. 2: Characteristic domain size after 87 ms of amplifi- 
cation at variable quadratic Zeeman shift g/. Data (circles) 
are averages over 5 experimental repetitions; error bars are 
statistical. Horizontal error bar reflects systematic uncer- 
tainty in g/. Predictions based on numerical simulations for 
[Aa| = 1.45 as '§J (squares) and 1.07 as [2l|] (triangles) are 
shown, with error bar reflecting systematic uncertainty in the 
atomic density. 

form condensate within the Bogoliubov approximation, 
one defines mode operators b n ,k — u4> n ,k + v<j)[ L _ k for the 
two polarizations of transverse (n € {x, y}) spin excita- 
tions. The spin-dependent many-body Hamiltonian H s 
is then approximated as representing two independent 
parametric amplifiers: 

H s = -t Y)E s {k)\{bl, k ~ C + h %M ~ • (2) 

k 

The parametric amplifiers serve to squeeze the ini- 
tial state in each spin excitation mode, amplifying one 
quadrature of b n j, (its real part) and de-amplifying the 
other. 

The above treatment may be recast in terms of spin 
fluctuations atop the paramagnetic state: fluctuations of 
the transverse spin, represented by the observables F x 
and F y , and fluctuations in the alignment of the spinor, 
represented by the components N yz and N xz of the spin 
quadrupole tensor. We identify the Bogoliubov operators 
defined above as linear combinations of these observables. 
Based on this identification, we draw two conclusions. 
First, an ideally prepared paramagnetic condensate is 
characterized by quantum fluctuations of the Bogoliubov 
modes. In the linear regime, fluctuations in b x ,k (b y ,k) 
correspond to projection noise for the conjugate observ- 
ables F x (F y ) and N yz (N xz ). Second, the dynamical 
instabilities lead to a coherent amplification of these ini- 
tial shot-noise fluctuations. While in the present work 
we observe only the magnetization, in future work both 
quadratures of the spin-mixing amplifier may be mea- 
sured using optical probes of the condensate nematicity 
[22| or by using quadratic Zeeman shifts to rotate the 
spin quadrature axes. 



FIG. 3: Temporal evolution of the transverse magnetization 
variance G(0), evaluated over the central 16 x 124 /im region 
of the condensate and averaging over 8 experimental repe- 
titions; error bars are statistical. The contribution to G(0) 
from imaging noise was subtracted from the data. Predictions 
from numerical calculations for |Aa| = 1.45 as and 1.07 as 
are shown as black and gray lines, respectively. 

To test the validity of this quantum amplification the- 
ory, we determine G(0) over the central region of the 
condensate after different intervals of amplification (Fig. 
[3]) . We consider the linear-amplification theory to be ap- 
plicable for t < 90 ms, and, following Ref. [11[, perform 
a least-squares fit to a function of the form 

G(0)| t = G(0)\ tm x y/tfae*-^ (3) 

Here r is the time constant characterizing the growth rate 
of the magnetization variance and G(0)| tm is the value of 
G(0) at time t m = 77 ms. 

To compare our measurements to the amplifier theory 
outlined above, we performed numerical calculations of 
G(0)| t , taking into account the the inhomogeneous den- 
sity profile; dipolar interactions; proper position-space, 
rather than momentum-space, spin excitation modes; 
and quantum fluctuations of the initial state 28]. Such 
calculations, results of which are shown in Figs. [3] and 2J 
were performed for several values of the scattering length 
difference Aa within the range of recent measurements 
[1,0] . Our data are consistent with the quantum- limited 
amplification of zero-point quantum fluctuations in the 
case that |A<z| lies in the upper range of its reported val- 
ues. Taking r to be determined instead solely on the 
basis of our measurements, the magnetization variance is 
measured to be between 1 and 50 times greater than that 
predicted by the zero-temperature quantum theory 

Several investigations were performed to identify tech- 
nical or thermal contributions on top of the expected 
quantum spin fluctuations of our samples. We place an 
upper bound on thermal noise by performing a Stcrn- 
Gerlach analysis of populations N± in the \m z = ±1) 
states just after the quench. Obtaining N± < 3 x 10 2 and 
assuming an incoherent admixture of Zeeman sublevels, 
the thermal contribution to G(0)|n is below (2N±/Nq) = 
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FIG. 4: The magnetization variance G(0)t,„ at t m = 77 ms 
and exponential time constant r of the amplifier, obtained by 
fitting data in Fig. [3] are indicated by contours of the 1, 2, 
and 3 a confidence regions using a \ 2 test. Predictions from 
numerical calculations of the zero-temperature quantum am- 
plification theory, assuming different values of |Aa|, are shown 
(circles and interpolating line). Error bars reflect systematic 
uncertainty in q/ and in the condensate density. Time con- 
stants corresponding to reported values for |Aa| are indicated 
at bottom. 



3 x 10~ 4 . Wc confirmed that our results are insensitive 
to variations in the gradient strength, duration, and ori- 
entation used during the initial state preparation, and 
also to the delay (varied between and 110 ms) between 
this preparation and the initiation of the spin amplifier. 
We checked against technical noise that would induce ex- 
trinsic Zeeman transitions during the experiment, find- 
ing that a condensate starting in the \m z = —1) state 
remained so for evolution times up to 400 ms following 
the quench. Altogether, these results suggest that para- 
magnetic samples were produced with a near-zero spin 
temperature, to which the quantum amplification theory 
may be expected to apply. 

The microwave fields used to vary the quadratic Zee- 
man shift cause a slight mixing between the F — 1 and 
F = 2 hyperfine levels. We checked for the influence of 
this admixture by varying the relative contributions of 
the static and modulated field contributions to the final 
quadratic Zeeman shift q/ = qs + q^- As expected given 
the small value of (fl/<5) 2 < 10~ 3 , no variation in the 
magnetization evolution for constant qf was observed. 

It remains uncertain whether the zero-temperature 
amplifier theory should remain accurate, out to a gain 
in the magnetization variance as high as 30 dB, in a 
non-zero temperature gas subject to constant heating 
and evaporation from the finite-depth optical trap. In- 
deed, previous work showed a strong influence of the 
non-condensed gas on spin dynamics in a two-component 
gaseous mixture 23[. We compared the amplification of 



obtained for different optical trap depths. We observed 
no variation, but note that the condensate fraction was 
not substantially varied in this comparison. 

In this work, we have demonstrated near-quantum- 
limited amplification of magnetization fluctuations us- 
ing dynamical instabilities in a spinor Bose gas. Just as 
the demonstration of matter-wave amplification in scalar 
condensates suggested a host of applications 13, IE 3] , 
the spin-mixing amplifier may serve as an important ex- 
perimental tool in probing systems with multiple degrees 
of freedom or in spatially resolved magnetometry [201 ] . 
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